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NOMENCLATURE 



A* throat area 

Ah throat area for High Speed CDTD 
Ai throat area for Low Speed CDTD 
a throat dimension 
b width test section 
D outer diameter 
e energy of the unit volume 
L length of blades 
m mass flow rate 
P power 

P3 specific power 
p pressure 
R radius of curvature 
ri inner radius 
r2 outer radius 
a blade spacing 
T torque 

U velocity component corresponding to 
Un velocity component normal to edge 
U2 wheel rim velocity 

V velocity component corresponding to 
Vt, velocity component tangent to edge 



coordinate 



coordinate 
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V 9i t.an 9 ential component of inlet velocity 
V ©2 tangential component of outlet velocity 
W2 relative exit velocity 
Wx work 

w' rotational velocity 
Z number of blades 
a efflux angle 

02 blade exit angle with respect to radial direction 
02' negative blade exit angle with respect to radial 
direction 

T specific heat ratio 
( internal energy 
"c power conversion efficiency 
^ density 

Notations used in Appendices B and C are independently 
defined therein. 
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INTRODUCTION 



I. 

The motivation for the work reported here came from two 
exploratory projects, a Rotary Detonative Wave Engine, 
(RDWE>, described in Ref. 1 and a high speed Centrifugal 
Diffuser Teat Device <CDTD), described in Ref. 2. Both 
projects required the geometrical design of an axisymmetric 
generator for nearly-tangential transonic outflow. 

To improve gas turbine engine performance an increase 
in cycle pressure ratio is needed. One alternative to 
mechanical compression is to supply the compression 
thermally by detonative combustion. If the pressure rise 
produced in this manner is sufficient for the engine, no 
compressor and associated drive turbine are required and a 
simpler more efficient engine might result. 

As described by Monks CRef. 13 , a proposed geometry 

for such an engine might include a rotor with outflow 
passages which are turned at the exit to produce nearly 
tangential flow. Figure 1 illustrates the proposed RDWE. 
Monks attempted to measure the torque on a static simulated 
rotor as a detonation wave propagated from the center 
outward. A reaction torque was produced from the nearly- 
tangential expulsion of gases from the rim. 
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The question then was, what would be the best design 
for the rotor passage to maximize the torque impulse? In an 
analysis included in Appendix A it was shown that tangential 
flow at the rim of the RDWE would yield the highest specific 
power and highest power conversion efficiency. A design was 
required therefore of nozzles which would direct the outflow 
almost tangentially. The nozzles would operate in an 
unsteady mode. 

The second project which motivated the present study 
involved the development of a test rig for radial 
diffusers. Such diffusers are used in centrifugal 
compressors to convert the kinetic energy of the flow from 
the rotor into static pressure. Erwin CRef. 21 proposed and 
designed a so-called low speed Centrifugal Diffuser Test 
Device, <CDTD) . The CDTD simulated the flow of air 
delivered to the diffuser from the compressor rotor using 
static swirling vanes. The facility would allow the 
evaluation of proposed new diffuser geometries and permit 
validation of computer analyses for diffusers operating in 
the fully subsonic flow regime. Figure 2 shows a schematic 
of the low speed CDTD. Vidos [Ref. 3] presented results 
from the flow generator used in the CDTD and performed an 
analysis of the internal flow. As mentioned by Vidos, the 
design of centrifugal diffusers is presently based on 
experimental results for two dimensional and conical 
diffusers. Numerical methods are currently under 
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development for viscous, three dimensional, transonic flow 
(a regime containing both supersonic and subsonic regions), 
with adverse static pressure gradients. Such methods are 
required to design optimum radial diffusers for transonic 
conditions. Equally, a high speed version of the CDTD was 
required in order to examine such diffuser flows under 
controlled conditions and to verify the emerging codes. 
Again, nearly tangential swirling outflow was required, this 
time in steady flow. 

The main interest here was therefore to design vane 
shapes and the apparatus for a high speed CDTD. The design 
would closely parallel the general arrangement of the low 
speed CDTD. The high speed CDTD, had the unusual 
requirement of radial-outflow with convergency at the rim 
exit. The limitation on the design was primarily that the 
hardware could be machined easily and inexpensively in- 
house. Computer programs for the design of supersonic 
nozzles for radial-inflow turbomachines were found to exist, 
but most were highly involved and would need extensive 
machining processes to generate the blade shapes. 

Particular designs are discussed herein. A literature 
search indicated that computer programs for radial -outflow 
devices using simple vanes for two-dimensional isentropic 
flow of a perfect gas were not available. 

In the present report, the design of the blading for a 
high speed CDTD model and the general arrangement of the 
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apparatus are described. General methods and computer 
programs for the design of double circular arc and wedge arc 
blading for radial -outflow vanes are described in detail in 
appendices. The report also discusses a method of 
predicting the flowfield through a selected passage using an 
Euler code based on the Godunov method of solution. 
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II. APPARATUS DESIGN 



Factors considered in the design of the high speed CDTD 
fell into two categories; namely, limitations due to the 
available air supply system and piping and the requirement 
that the test section have dimensions appropriately scaled 
for the diffusers required to be tested. 

A. DESIGN CONSIDERATIONS 

The high speed CDTD was designed to operate as an 
intermittent blow-down wind tunnel. In order to use the air 
supply system available in building 230 at the Naval 
Postgraduate School, the test section flow area was dictated 
by the available mass of stored air and the minimum desired 
runtime of two minutes. Using proportions similar to those 
used in the low speed CDTD (see Figure 2), the radius of the 
swirl generator became four inches while the radius at the 
test section was five inches. The method used for 
calculating these dimensions is given in the following 
section. The difference between the swirl generator outer 
radius and the test section radius allowed for mixing of the 
wakes from the blading (or nozzles) used in the swirl 
generator. The swirl generator was limited to a length of 
two inclies to meet the scaling requirements set by the test 
section. Because of the larger pressure differential 
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required across the swirl generator in the high speed 
version, turbine like turning vanes were selected rather 
than the sheet metal nozzle construction used in the low 
speed device. 

B. APPARATUS DESCRIPTION 

Figure 3 illustrates the overall arrangement of the 
proposed CDTD design. The swirl vanes for the high speed 
CDTD are sandwiched between two circular plexiglas disks. 
Plexiglas disks were used to allow visualization of the flow 
through the vanes. The disks are mounted on an aluminum 
center pipe shown in Figure 3. The pipe provides double 
axial in-flow. Following the work performed by Feieresien 
CRef. 4] a transition section to turn the air from the axial 
to the radial direction was added in the form of a plate on 
a tie-rod on the pipe centerline. The design of the 
transition section is not complex if the area of the radial 
outlet is smaller than that of the total axial inlet. An 
acceleration through the transition section will then occur, 
leading to generally favorable pressure gradients. These 
consideration required a double-entry supply arrangement for 
the CDTD model. 

The piping is attached to the air supply system. Figure 
4 shows the configuration of the blades on the plexiglas 
disks and in relationship to the supply pipe. 
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C. DESIGN CALCULATIONS 

Vavra CRef. 53 showed that in the design of a turbine 
cascade the efflux angle a can be calculated from the 
geometry of the cascade and that it is independent of the 
incidence angle i of the inlet velocity, see Figure 5. 

Vavra also noted that it is common practice in steam turbine 
design to use the relation 

a = arc cos <a/s) (1) 

where •*a” is the throat opening at the discharge and s is 
the blade spacing. Equation (1> applies when the trailing 
edges of the blades are of negligible thickness. In 
practice the "wedge** angle at the trailing edge and trailing 
edge bluntness modifies this dependence. 

With the limitation of a two minute runtime and the 
available mass of stored air. Demo [Ref. 63 found that for 
the supersonic cascade rig in building 230 an area of 6.27 
square inches was needed at the nozzle throat where the flow 
would choke. This same area was taken for the total outlet 
throat area of the swirl vanes for the high speed CDTD. 

Referring to the following sketch, the area at the exit 
where the flow will choke, A», is given by, 

A« = L a 2 (2) 

where L is the vane length, a is the throat width and Z is 
the number of vanes . 
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Figure II. 1 Description of Conditions at Throat 



The blade spacing, s, is given by 

s » n D / Z <3) 

where D is the diameter at the vane exit. Froa equations 
(2> and <3) , 



A- ■ L (a/a) n D. (4) 

In Ref. 7, the following conditions were calculated for 
a “full scale*' high speed CDTD design to have a test section 
width of 1 inch at 2S inches diameter: 

at Qi»55o, at Ri*12.5 inches 
Vial392 ft/aec at Tq* 520 «R, giving 
Vei»1140 ft/sec and Vro» 1357 ft/sec and for 
a swirl vane length of 5 inches, VRo*‘i®2 ft/sec 
and Bo”d2.34«. An increase of swirl vane length 
to 5.713 inches gave Vro« 160 ft/sec and Ro»03.3« 



Using the above data for the “full scale" high speed 
CDTD in equation (4), A* « 3d. 5 inches^. Using a linear 
scaling to the model in which A» ■ 6.27 inches^, the scale 
factor is 
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a (6.27/38.5) S: .40. (5) 

Uelng this scsle fsctor, ths dimensions cslculoted for the 
high speed model sre shown In the following sketch. 




Figure II. 2 Cross Sectional View of the Test Section 
From equations (2)> 

a z » A* / L (6) 

therefore for a diameter of 8 Inches we can predict a 
required throat dimension, and blade spacing using equation 
(3) for a selected number of blades. Results are given In 
table 1. The requirement then Is to design suitable turning 
vanes to produce the required throat area and spacing. 
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Ill, BLADING DESIGN 



A. BACKGROUND 

In Ref. 5 Vavra gave a procedure to establish a blade 
profile for specified values of throat width a, spacing s, 
trailing edge thickness te# and inlet flow angle, a. Figure 
6 illustrates these features. 

Vavra^s design was for axial turbine blades or vanes. A 
computer program was available to construct blading 
following Vavra's methods CRef • 5] • The method used 
circular arcs, and straight lines. Wedge angles and leading 
and trailing edge radii were input parameters to determine 
a flow channel with suitable acceleration. No comparable 
method and program were found for radial outflow vanes. 

In 1968 Goldman and Scullin CRef. SI designed blades 
with inlet transition arcs, circular arcs and outlet 
transition arcs. However, these blades were not converging 
at the exit and were not designed for radial-outflow. 

Later Goldman and Vanco CRef. 93 developed a computer 
program for the design of sharped-edged throat supersonic 
nozzles for isentropic flow. Their program was constructed 
for the radial-inflow case. The Goldman and Vanco design 
also involved the requirement of contour machining. 

There is a substantial amount of information on both 
radial-inflow and axial turbine blade design in the 
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literature CRef. 8 , 3 , 10 , 11 , and 12 for example]. Yet none 
was found which dealt with the required case of radial- 
outflow, 

B, DESIGN METHOD 

The expansion process in a radial turbine differs 
appreciably from that in an axial turbine because of the 
radius changes, Vavra^s work CRef, 5] involved axial 
turbine design. The work presented here for radial-outflow 
vane design followed an appreciation of Vavra's methods and 
possibly retraced steps of others who have previously sought 
the simplest geometries for machining purposes. 

The simplest shape^ that of double circular arcs was 
attempted first. Its development is presented in Appendix 
B, It was anticipated that with the requirements of radial 
inlet flow and tangential exit outflow the double circular 
arc method could be used in the design of radial-outflow 
turning vanes. Once the method was developed and exercised 
it was clear that more geometrical constraints were required 
to produce a nozzle with an exit which was both nearly- 
tangential and convergent. 

The double circular arc program generated vanes which 
appeared to converge within the passage and then diverge 
towards the exit of the nozzle. Therefore, the flow would 
choke within the nozzle and the outlet flow conditions could 
not be predicted very well. 
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Modifications were made to the double circular arc 



design to achieve convergency at the exit. The first step 
was to remove the pressure line and use the camber line as 
the pressure surface. The assumption was made that the 
effects associated with inlet flow incidence resulting from 
the removal of the pressure line would be negligible. The 
results indicated that the passage was indeed converging 
toward the exit. The second step in an attempt to assure 
tangency and convergency at the exit was to allow for a 
thickness and a straight section at the leading and trailing 
edges. The straight section was such that a specified 
throat dimension was achieved at the trailing edge of the 
adjacent previous vane. The resulting program^ which is 
presented in Appendix C, designs a wedge shaped arc blade 
such that a passage which converges at the exit is realized 
with a specified throat dimension. Using the method a high 
speed CDTD blading design was obtained. Figure 7 shows the 
resulting blade design. 
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IV. FLOW ANALYSIS 



A. METHOD 

A FORTRAN computer program for solving numerically the 
two-dimensional Euler equations using the Godunov method was 
developed at the Turbopropulaion Laboratory by Dr. S. 
Eidelman [Ref. 121. With Dr. Eidelman^s help the program 
was to be applied to obtain a prediction of the flowfield 
through the vane passage. An abbreviated description will 
be given here of the computational method. 

The unsteady two-dimensional Euler equations are written 
as 

^ t hx by 

where 







pu 




pu 


u 


, b = 


p + P U2 


# C = 


p UV 


V 




UV 




p+ P V2 


e 




<e+p)U 




<e*p) V 



2 2 

e = ('<('►(11 '►V)/2) = energy of a unit of volume 

( = P = internal energy 

< T -1 )|0 

^ = density 

U and V are the velocity components corresponding to the x 
and y coordinates 
p = pressure 
T = specific heat ratio. 
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Equation (7) is nondimensional ized using arbitrary 
reference quantities, and the Godunov method [Ref. 123 is 
applied to compute at each time interval the fluxes at the 
boundaries of the cells of a grid established for the 
problem under consideration. The steady state solution is 
found when no further change occurs between successive time 
steps . 

The computational domain is shown for a typical cascade 
problem in Figure 8 and the grid and cell notation are 
defined in Figure 9. Referring to Figure 8, for the 
tangency condition at the wall of the passage a solid wall 
boundary condition is imposed over the segments 3-5 and 4- 
5. Boundary conditions for segments 1-3, 5-7, 2-4 and 8-8 
are for periodicity of the flow with respect to the y axis. 

The parameters held constant at the inlet are the x 
components of velocity, U, the flow angle, 0, the gas 
stagnation enthalpy, H, and the entropy, S. The velocity, 
density and pressure can be calculated from these quantities 
easily as ^±xip Pin# ^nd Pin* At the outlet boundary only 
pressure is defined and other flow parameters are calculated 
in the solution. In the solution procedure, Uout# ^out# 

P out s®t equal to the values of U, V and p one point 

ahead of the outlet boundary. 

A subroutine (program) constructs a non-orthogonal grid 
covering the computational domain as shown in Figure 9. The 
flow parameters related to the center of the cell have a 
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fractional index: 3 + h/ The cell boundaries have one 

fractional and one integer index: j, or The 

parameters at the time t will have subscript indices and at 
time t-^ St the parameters will have superscript indices. 

The Euler equation is approximated to first order by the 
equation in table 2. Solution to the Riemann problem is 
initiated assuming that the left state is that located at 
and the right state is that at 3+)^,k+h. 

The solution to the Riemann problem will give 

Un • the velocity component normal to 
the edge < 3 ,k*^l) -< 3 #k) 

Vt = the velocity component tangent to 
the edge ( 3 ,k+l) -< 3 /k) 

The capital letters R (rho)^ U, P are the parameters 

calculated at the edges of the zones. 

The component of the velocity is then transformed back 
to the cartesian coordinate system and a final value of the 
parameter on the cell edge is received 

<R^P^Un/B) , , 

3 h 



B. RESULTS 

The program (EGE) was run on the IBM 370-3033 at the 
Naval Postgraduate School. Figure 10 illustrates the grid 
covering the passage for the wedge-arc blading shown in 
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Figure 7. An enlargement of the exit area is shown in 
Figure 11. 

In order to apply the code to the radial geometry, some 
modifications were necessary- A transformation of 
coordinates was made to conveniently accommodate the radial 
flow inlet boundary condition. Since the program remained 
basically two-dimensional, the outlet flow was treated as 
being periodic over the blade space at the exit which then 
did not change with radius. Thus the solution found will be 
valid for the flow within the rim and approximately true 
just outside the rim. For this reason the computational 
domain in Figure 11 was not extended radially outwards 
beyond one blade space. 

Results could not be obtained for the flow through the 
passage within the time available. Some dif f iculities were 
encountered by the program in the region of the throat. 

Some modifications of the grid must be made before a 
converged, steady flow solution can be obtained. 
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V 



SUMMARY AND DISCUSSION 



Two exploratory projects at the Turbopropulsion 
Laboratory^ 1). to teat the concept of detonative combustion 
aa an alternative to mechanical compression in gas turbines 
and 2) . to design a teat rig for radial diffusers used in 
centrifugal compressors, were examined. The projects were 
similar in nature since they involved the generation of 
nearly-tangential radial outflow and each required a design of 
turning vanes for a test apparatus. The overriding limitation 
for either project was that any apparatus developed would have 
to be easily and inexpensively manufactured in-house. 

In the Rotary Detonative Wave Engine, "RDWE”, project, 
turning vanes were required such that upon detonation the 
maximum possible torque impulse would be obtained. The best 
design possible to yield maximum power hence maximum torque 
was to produce turning vanes that exit nearly- tangential . 

Using Euler's Turbine equation a relationship was developed 
relating exit angle to specific power and efficiency. 

The results for the unsteady expansion process showed 
that the value of the relative flow angle (B2'> established the 
limiting conversion efficiency and specific power. If the 
relative velocity increases beyond the limiting point the 
specific power generated will be absorbed rather than 
generated and efficiency will decrease. High efficiency 
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occurs in large part only when the specific power is low or at 
the limiting point. In general the conversion efficiency will 
vary during the expansion process, but will always be poor. 

This result implies that the concept of the RDWE can not 
be efficient as a stand alone rotor. Other methods would have 
to be employed to increase efficiency, such as the utilization 
of additional components, particularly those which would allow 
nearly constant flow by sequencing the detonation wave with 
respect to partial ports. 

The second project involved the development of a test rig 
for a high speed Centrifugal Diffuser Test Device, **CDTD**. 

The project would closely parallel an in-house concept of a 
low speed CDTD. The low speed version simulates the flow of 
air delivered to a diffuser as a steady process. The high 
speed version simulates in a steady process, the transonic 
Mach numbers characteristic of advanced centrifugal 
compressors . 

For the high speed CDTD a passage would be required which 
would cause the flow to converge and choke at the exit of the 
passage. Because the concept of the RDWE was questionable the 
main interest tuned to designing vane shapes and the apparatus 
for a high speed CDTD. 

Two methods and their respective computer programs were 
developed. The first method involved the development of 
double circular arc vane shapes. Given the outer radius, the 
inner radius and the number of blades the radii of the suction 
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and pressure sides were determined. The design emphasis was 
on acheiving a tangent suction surface at the exit and a 
radial camber line at the inner radius. The concurrent 
condition of a tangent camber line at the exit gave the 
solution of a blade of zero thickness. 

The geometrical description resulting in a tangent camber 
line at the exit, resulted in the sine of the camber angle 0 
being equal to the cosine of the subtended blade angle . 
Following manipulations to relate the value of ^ to just the 
radius ratio (the only specified parameter), an equation of 
fourth order in sin resulted. Using Newton's method of 
interation the first (smaller) root was found. Using this 
value of a blade shape was generated. It was clear that 
the suction aide was tangent but the efflux angle a was not 
zero. An efflux angle of zero would indicate that the camber 
line was also tangent to rim. A graph of the fourth order 
function indicated that two roots were possible for each 
radius ratio. Only the first root was solved for. It is 
probable that the second root would yield the tangency 
condition of both camber line and suction side. The selected 
root was shown however to give sin 0 - cos 'f" . 

It was realized that convergence at the exit could not be 
guaranteed with the double circular arc blade design. 
Modifications were applied to the design to cause the passage 
to converge. Improved convergence resulted when the original 
camber line and the suction line were used to describe the 
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blade 



The simple double circular arc design required curved 



or thin trailing edges. With the addition of a finite 
thickness at the leading and trailing edge, a straight line 
segment, and a wedge angle all of which would better direct 
the flow in the tangent direction, practical blade shapes with 
convergent passage were acheived. 

Geometrical equations were developed to yield the proper 
blade design. The exit throat dimension was fixed in the 
design based on the limitations of the available air supply 
and based on dimensions scaled from a high speed CDTD design. 
It was anticipated that the flow would choke at the specified 
throat dimension between adjacent auction and pressure sides 
and the straight line segment would intersect the exit circle 
which was tangent to the suction side. This did not occur. 

The length of the straight line segment varied between 1 and 5 
percent of the value needed for this condition. It was not 
required that the exit circle be tangent to the pressure side. 

From the geometry and the specified parameters a blade 
was developed giving convergency at the exit. The camber 
angle 0 and the subtended blade angle ^ were specified from 
experience obtained in the development of the double circular 
arc design. The chosen angles 0 and ^ gave an efflux angle a 
between 8 and 9 degrees. To achieve the perfect condition of 
tangential flow angle a would be required to be equal to 
zero. Based on current technology and manufacturing abilities 
an angle in the area of 8 or 9 degrees would give as close to 
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tangential flow as is possible. It is noted that the angle 0 
is related uniquely to the angle . Once the radius ratio is 
given • 

The coordinates of the designed passage were then input 
into the FORTRAN program (EGE) to predict analytically the 
flowfield through the vane passage, by solving the Euler 
equations. The program (EGE) was designed to run on the IBM 
370-3033. The program successfully constructed a non- 
orthogonal grid covering the domain, however results could not 
be obtained for the flowfield through the passage in the 
initial attempts. It 

is thought that the difficulty lies in the curvature and 
density of the grid in the region of the exit throat. Since 
the EGE program has not previously been applied to a radial 
outflow problem with severe curvatures the task of analyzing 
the present flow with the EGE code is now thought to require 
the concentrated effort of the code's author. 

In contrast, the CDTD vane design can be manufactured 
without difficulty and the flow from the turning vanes can be 
measured experimentally. It is suggested that such an 
experiment would provide a severe test of the ability of the 
EGE code to simulate transonic radial outflow containing 
embedded waves. 
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VI . CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The following conclusions are offered: 

The Rotary Detonation Wave Energy concept in its original 
configuration is potentially inefficient. An increase in 
efficiency would require the addition of a component which 
could recover the unsteady kinetic energy of the flow leaving 
the rotor. 

A program for designing radial outflow vanes using double 
circular arc blades was demonstrated successfully. While the 
method may be useful in other contexts, it was found to be 
unable to produce converging passage shapes with a near- 
tangential outlet flow direction. 

A program for designing radial outflow vanes using wedge- 
arc blade shapes was used successfully to produce converging 
passage shapes with prescribed throat widths and with a near- 
tangential outlet flow direction. 

Both blade shapes can be produced relatively inexpensively 
from standard metal stock. 

A preliminary design of a test model of a high speed 
Centrifugal Diffuser Test Device was produced. 
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B. RECOMMENDATIONS 



The following recommendations and areas of further 
investigations are suggested: 

An arrangement of partial porting and following power 
turbine for the Rotary Detonation Wave Engine should be 
examined • 

The program for double circular arc blades should be 
modiied to seek both roots of the condition for camber line 
tangency, and to offer the selection of any value of between 
them. This will result in the ability to examine any value of 
9^ within the range of any practical interest. 

In the case of wedge-arc blades^ the possibility of 
solving for the particular case in which the circle at the 
passage outlet is also tangent to the pressure side should be 
examined. A procedure should also be introduced to calculate 
the progression of the passage area taking into account the 
increasing radius . 

The proposed blading geometry should be constructed and 
the outlet flowfield examined exper imently . 

Work should be continued to obtain a computational 
prediction of the flow through the passage. The model test 
could then serve as a verification of the code. 
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TABLE 1. Required Deaign Dimensions for Convergent 
Throat in High Speed CDTD . 



NUMBER OF BLADES 
Z 


THROAT DIMENSION 
a < in . ) 


SPACING 
s ( in. ) 


6 


.5225 


4.1888 


12 


.2613 


2.0944 


18 


.1742 


1.3963 


36 


.08708 


.6981 


45 


. 06967 


.5585 


60 


. 05225 


.4189 
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TABLE 2. First Order Solution of 
Euler Equations 



- j+h,k+»5 _ - 



= a 



j+h,k+Js 



At 



j+h/k+h 



' ^^j+h,k ” ^j+h,k ^j+h,k 



®j+>s,k+h ^^j + l,k+h “ *“j+l,k+Js * '^^j + l,k+is 



^j+h#k+l * ^^j+h»k+l ” ^j+*5/k+l '^^j+h/k+l 



^j,k+h ^j,k+h ” ^j/k+h ^^j,k+h^ 



where 

'^j+^/k+h ~ ^ ^^j + l,k+l ” ^j/k^^^j,k+l ” ^j + l,k 

(^j,k+l - ^j+l,k)<^j+l,k+l - Yj,k) 



^jth,k = ^j.l,k-^3,k 





p 




RU 




RU 


a = 


pu 


, B = 


P+RU^ 


, c » 


RUV 




pv 




RUV 




P+RV^ 




e 




(E+P) U 




(E+P) V 
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Flgur* 1. Schematic of Datonatlon Rotor Concept. 



Sv'lrl Vane Cylinder 
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Figur# 2. Schanatic of tha Low Spaad CDTD 
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Figure 3. SchesiQbic o£ the Proposed High Speed CDTD. 




Note: Not to scale. 



Figure 4. CDTD Flow Generator - Section AA 

from Figure 3. 
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Flgur« 5. Two-Dimensional Turbine Cascade. 
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Figur* 6. Axial Turbine Blade Geoeetrical Design. 
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Ro *4 Ri • 3.2 No. Bladea ■ 60 

9.1.72521902 0 ■ 80.82738092 

Condition of Tangency ia 
ain 0 = coe 

sin 0 = 0.159407735 cos » 0.159409430 
chord ■ 0.983540343 chord test ■ 0.983540073 
alpha =» 18.34514098 
radial solidity ® 1.909859317 
thickness to chord ratio ■ 0.194933797 
radius of camber circle ■ 0.758551242 
radius of suction circle ■ 0.575175292 
radiua of pressure circle ■ 1.462802226 



Figure 7. Final Blading Design. 
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Flgur« 8. Computational Grid and Boundary of the 
Computational Domain. 
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a) Grid Location in Coordinate Syate*. 




b) Single Cell 



Figure 9. Coeputational Grid Notation. 
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Fl9urtt 10. Computational Grid for Outflow 

Analyai« 
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Figur# 11 • Enlargement of the Grid in the 
Throat Region. 
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APPENDIX A 



TORQUE PRODUCED BY DISTRIBUTED GAS REACTION 
A.l METHOD 

The RDWE concept involves an unsteady process. The 
potential performance of the RDWE can be understood , 
approximately, by examining the range of conditions through 
which the unsteady process passes. Thus the steady flow 
performance of the proposed arrangement needs to be 
examined. We first attempt to predict analytically the 
power or torque impulse that can be generated from the 
engine under steady flow outlet conditions. 

The passages in the RDWE exit nearly tangentially 
(Figure 1). Since this condition is difficult to achieve we 
examine the effect of the exit angle on the amount of power 
that can be generated and the efficiency that can be 
produced . 

The basic equation for relating the fluid flow 
conditions to the energy transferred to the shaft is Euler's 
Turbine equation. The torque of the fluid acting on the 
rotor is given 

T = m < ri V01 - r 2 V02) Ad) 

If the rotor has a rotational speed w', the shaft power 
transferred or rate of work transferred from the fluid to 
the rotor is given by. 



49 



p 



w" T 



A<2) 



The work per unit mass flow is given by 

P / m = w' (ri Vei - r2 Ve2> A<3) 

The subscript 1 denotes entrance station, 2 denotes exit. 
Defining the rim speed as U = rw', equation A<3) becomes 

P / m = Wx = Ui Vei -U2 V02 A(4) 

A general velocity diagram for a rotor is shown in Figure 
Al. 

The tangential component of velocity VO is taken 
positive in the direction of the rim velocity U. Figure Al 
shows all positive velocity components in conventional 
cylindrical coordinates . 

The gas enters the wheel near the center, effectively 
with zero moment of momentum (about the axial direction) • 
Therefore V01 = O, and from equation A<4), 

Wx = P = -U2 V02 A<5) 

n 

For Wx to be positive, it is seen that V02 is negative. The 

exit velocity diagram representing the flow at the rim of 

the wheel is given in Figure A2 . 

From the geometry shown in Figure Al, in general 

sin ft = V0 ■ U A<6) 

W 
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Figure Al. Velocity Triangle of Conventional 
Cylindrical Flow. 





Figure A2. 



Velocity Diagram at Exit of Rim. 
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therefore 



V© = W sin B + U A<7) 

Substituting equations A(7) into A<5) we obtain, 

Wx = P = - U 2 C W sin B ♦ U ] 
m 

= - U^2 C_W_sin B ♦ 1 3 

U2 

= U^2 C W 3 in <-B> - 1 ] A (8) 

U2 

Defining the power conversion efficiency ^c as the 
fraction of the available power actually converted to shaft 
work, thus 

= P/m A<9) 

p/m + v2^/2 

and defining B2'= “*^2^ relative flow angle at the rim 

which is positive in the direction of rotation, then 

P = U2 C W 2 sin B 2 ' - 1 3 . A<10) 

m U2 

From the geometry shown in Figure Al, 

V2^ = V©2^* < W2^ -CV©2-U23^) 

SO that, 

W2^= CU2 V©2 - W2^3 A<11) 

Using equations A<5> and A<10) 

- U2 V ©2 = -Wx = C_W 2 sin B 2 ' - 13 A<12) 

U2 

Substituting equation A(12) into A(ll) gives 
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= h U 2 ^ [ - 2 W? sin O 2 ' ♦ 1 ♦ W-? ^ ] A<13) 

U 2 U 2 

Substituting equation A<10) and A<13) into A<9) 

"c = P/m 

P/m hV2^ 



' U_ (W_/U_) sin - 13 

_2 — 2 — 2 2 r 

U^[<W 2 /U 2 >sin B 2 ' - 13 + !tU 2 ^ £-2Wsin(32' *1* W 2 ^3 



U 2 



U 2 



= (W-?/Uo) sin - 1 

-a * h <W2^/U2^> 



»c = - 2 C <W-?/U? sin B 2^ - 1 3 
< 1 - W2/u2) 



A(14) 



Equations A(10) and A(14) are valid for a discreet jet, 
or for distributed gas at the rim. For the distributed gas, 
it is possible to calculate the referred power as a function 
of wheel radius, giving specific power, P 3 . Using equation 
A<10) , 

Pa = P/mU2 = UW 2 /U 2 ) 3in B2'- 13 A<15) 

Therefore equation A (14) becomes, 

"<= = 2 Pa 

<W2/U2>2 - 1 

For sin B 2 '= 1» implying that B2'=90 degrees, 

t>c = 2 

(W2/U2) ♦! 
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Thus the conversion efficiency is unity when W 2 /U 2 = 1 
but this is a condition which, from equation A(15) gives 
zero power output. 

A. 2 RESULTS 

The value of the specific power and conversion 
efficiency (given by equation A(15) and A(16)) as a function 
of the ratio of relative flow to wheel velocity, for various 
values of the relative exit angle, are shown in Figure A3 
and Figure A4 respectively. 

Figure A3 shows that the specific power increases 
linearly as the exit velocity increases, and decreases as 
the relative exit angle becomes less tangential . The 
conversion efficiency is seen to decrease as the power 
increases and to decrease as the relative flow angle 
decreases (less tangential) • 

The results shown in Figure A3 and Figure A4 can be 
used to gauge the performance of the unsteady RDWE. In the 
unsteady expansion process the value of the relative flow 
angle (fl2^ will be constant, being determined by the passage 
exit geometry. The wheel speed will be approximately 
constant and only the value of the relative velocity will 
change. Thus the unsteady process at different times will 
be described by points on a single line shown on Figure A3 
and A4, depending on the particular value of R 2 • 
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VELOCITY RATIO VS SPECIRC POWER 




y3MOd Dldl03dS 
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VELOCrrY RATIO vs CONVERSION EFRCIENCY 




AON3IOUd3 NOISU3ANOO 



56 
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RELATIVE VELOCnr / RIM VELOCfTY 



The difficulty, clearly, is that the exit velocity must 
not be allowed to decrease below a specific factor times the 
wheel speed or power will be absorbed rather than generated. 
Also the conversion efficiency will vary during the 
expansion, but will always be poor. 

This implies that the concept of the RDWE as shown in 
Figure 1 can not be efficient. The device would require 
some type of diffuser that would convert the kinetic energy 
leaving the rotor into static pressure. Unfortunately, a 
variation in the relative velocity magnitude implies 
variation in both magnitude and direction of the exit 
absolute velocity. A suitable diffuser for such a flow does 
not exist. Consequently, other geometries, particularly 
those which would allow constant flow to the diffuser by 
sequencing the detonation wave with (partial) ports at the 
rotor exit, should examined. 
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TABLE A1 



Specific Power and Power Conversion Efficiency 
Va Velocity Ratio Data 





60 DEG 


"C 


70 DEG 
Wo/U^ P« 




0 


-1 


2 


0 


-1 


2 


.5 


- . 56699 


1.51200 


.5 


-.53015 


1.4137 


1.0 


- . 13397 


0} 


1.0 


-.06031 


CD 


1 .5 


. 29904 


.47846 


1.5 


. 40953 


. 65526 


2.0 


. 73205 


. 48803 


2.0 


.87938 


. 58626 


2.5 


1.16506 


.44383 


2.5 


1 . 34923 


.51399 


3.5 


2.30109 


.36108 


3.5 


2.28893 


.40692 


4.5 


2.89712 


.30099 


4.5 


3.22862 


.33544 


5.5 


3.76314 


.25731 


5.5 


4.16831 


.28501 


6.5 


4.62917 


.22444 


6.5 


5.10800 


.24766 


7.5 


5.49519 


. 19892 


7.5 


6.04770 


.21892 




80 DEG 




90 DEG 




W^/Uo 


P« 


nc 








0 


-1 


2 


0 


- 1 


2 


.5 


- . 50760 


1 . 35359 


.5 


-.5 


1 . 33333 


1.0 


-.01519 


OD 


1.0 


0 


1 


1.5 


.47721 


.76354 


1.5 


.5 


. 80000 


2.0 


.96961 


.64641 


2.0 


1.0 


.66667 


2.5 


1 . 46202 


. 55695 


2.5 


1.5 


.57143 


3.5 


2.44683 


. 43499 


3.5 


2.5 


.44445 


4.5 


3.43164 


. 35653 


4.5 


3.5 


. 36364 


5.5 


4.41644 


.30198 


5.5 


4.5 


. 30769 


6.5 


5.40125 


.26188 


6.5 


5.5 


. 26667 


7.5 


6.38606 


.23117 


7.5 


6.5 


.23529 
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APPENDIX B 



DOUBLE CIRCULAR ARC BLADING DESIGN 



B,1 METHOD 

A computer procedure was developed to design radial 
out-flow, near tangential turning vanes, specifying only the 
outer radius, the inner radius, and the number of blades. 

The blade spacing, radial solidity, camber angle, suction 
and pressure surface radii, chord length, thickness to chord 
ratio, and efflux angle are determined. The suction surface 
is designed to be tangential at the rim. 

A Hewlett Packard 9830A computer was used. The blade 
shapes were generated using a Hewlett Packard 9862A plotter. 

The geometrical approach to the design of the vanes is 
illustrated in Figure B1 and much of the notation is defined 
in this figure. The following assumptions were made as 
initial constraints: 

1. radius at exit Rq^ is known 

2. radius at entrance Ri, is known 

3. number of blades 2, is known 

4. The suction side is tangent to the outer 
radius (point C) . 

5. The camber line c^, is tangential at the inlet 
radius (point A) . 

6. The blade subtends an angle at the center. 
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Initially, a procedure was developed which required 
specifying the camber angle 0, Rq, in addition to the 

number of blades. Since the choice of 0 which would give 
acceptable values of efflux angle and blade thickness was 
not obvious, an equation for the condition that the camber 
line was tangent at the rim (which would give a blade of 
zero thickness, since the suction surface was also required 
to be tangent) was written. The equation was seen to have 
two possible roots. The selection of the smaller 
(incorrect) root gave positive values of the efflux angle 
and reasonable values of blade thickness. The program which 
was developed provides the option of selecting 0, or of 
selecting the blade which results from finding the incorrect 
smaller root to the tangent camber line condition. 

Described here are the equations for the generation of the 
blade shapes. 

Defining peripheral blade spacing s as, 

s = 2 n Ro / Z 

and the angle subtended by the blade apace as 

/s = 360 / Z 

the radial chord as 

^ = (Ro “ Ri> / Ro 

then the radial solidity is given by 

(Tr = <Ro - Ri> / s = p 
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Using t-hese expressions and Figure B1 the following 



equations can be written! 

GE ~ ^c ~ ^c cos ( 0/2) 

= Rc Cl “ cos ( 0/2) 

HE = c T/2 ♦ Rc C 1 - cos < 0/2] 

Where T is the thickness to chord ratio. Also 

HE = Rs C 1 - cos (0S/2)] 

To solve for the thickness and T we subtract the 
equations B(l> and B(2) 

HE - GE = t/2 = T c/2 

Using Figure Bl, 

FE = Rc Cl - cos (0/2)3 - T c/2 

= Rp Cl - cos (0p/2)3 

Subtracting equation B<5) from B(l), 

GE - FE = t/2 = T c/2 

The chord c, is given by, 

c = 2 Rp sin (0p/2) 

= 2 Rs sin (0S/2) 

= 2 Rc sin ( 0 /2) 

All equations were made dimensionless by dividing by 
giving for example 

C/Rq = c" 

c“ = 2 Rp. sin ( 0p/2) 

Ro 

= 2 Ro s in <0s/2) 

Ro 



B(l) 

B<2) 

B<3) 

B<4) 

B<5) 

B<6) 

B<7) 

B<8) 

B(9) 

B(IO) 

o » 

B<11) 

B(12) 

B(13) 
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B(14) 



= 2 sin < 0 /2) 

For the geometry we seek, given Rq# and the camber 
angle 0 and efflux angle a, the angle AGO <or A) is given by 

A = 90 - < a ♦ 0/2) B(15) 

From the triangle AGO 

}>' * so - < a 0/2) 180 - 0/2 = 180 

resulting In, 

SO f = 0 ♦ a B(16) 

This equation fully defines a. The values of ''f' and 0 
are needed to complete the description of the camber line. 

To solve for the radius of each blade arc the following 
relations were used. 

First with 



As = Ec ” Ec cos 0/2 = HE 
Ro ^o 



B(17) 



and 



Ap - Rp “ Rp cos < 0p/2) = FE B<18) 

Ro Ro 

Using equation B<3), B<6) and B<7) 

Rc C 1 - cos <0/2)3 - Rp C 1-cos <0p/2)3 = T c/2 B<19) 

Also, 

As - T c/2 - Rp Cl - cos <0p/2)3 = T c/2 B<20) 
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which gives 



Ap = As - T c B<21) 

Using this relationship between As and Ap we can solve for 
Rs and Rp. From equation B<2), and B<3) , 

As = Rs - cos <06/2) B<22) 

or 

cos < 0s/2) * 1 - As/Rs 

Squaring and adding together equations B(22), and B<10> 
we obtain 

COS^ (0S/2) ♦ Sin^ < 0S/2> = <1 - Ag/Rg)^ + <C/2Rg)^ 
or 

1 = 1- 2<As/Rs> ♦ <As/Rs>2 -k c2/4Rs 2 
which becomes 

Rs = HAs ♦ c ^_1 B<23) 

8As 

or 

Ra = J4 As c ^ 1 B<24) 

Ro Rq 8As 

Using equations B<5), B(6), and B<16) a similar manipulation 
will yield, 

Rp = h Ap > c^ 1 B<25) 

Ro SAp 

The angles 0s and 0p are found using equations B<13), 
and B<14) since 

sin <0s/2) = c/2Rs 
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or 



0s = 2 arc sin ( c/Ro) 

2R3/R0 



B<28) 



and similarly 



sin <0p/2) = c/2Rp 



or 



0p = 2 arc sin ( c/Ro ) 

2Rp / R q 



B<27) 



The desired condition is that the suction surface be 
tangent to the outer rim. From Figure B1 this is given by 
the relation 



r is not an independent parameter since it can be 
calculated using equation B<25), thus 

r = C <Rcx C 1 - cos ( 0 s/2)3 - Ri^ [ 1 - cos ( 0 /2)3 ) B<29) 

(c/2) 

To insure near tangential flow at the exit we first 
seek the condition of a tangent camber line at the exit. 

This implies since we also seek a tangent suction surface, 
that the blade would have zero thickness. Using the 
notation in Figure B2 the following equations apply: 



0S/2 = a 0/2 



B(28) 



sin ( 0/2) = L / c 



BOO) 



and 



3in = L/ Ro 



B<31) 



Substituting equation BOO) into BOl) gives. 



c = Rr. sin ^ 



B<32) 



sin <0/2) 
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From the geometry 



cos ( s5/2) = y/c B<33) 

and 

cos f = R< ♦ V B<34) 

Ro 



or 



Ro cos = Ri y 

Substituting equation B(33) into B<35) gives, 
Ro cos = Ri + c cos < 0/2) 



B<35) 



B(36) 



or 

Ro cos "P - Ri = c cos (0/2) B(37) 

Equation B<32) can be rewritten such that, 

Ro sin P - c. sin (0/2) B(38) 

and dividing equation B(38) by B(39), 

R,^ sin ^ = tan (0/2) B(39) 

Ro cos ^ - Ri 

This equation establishes the relationship between^ 
and 0, Using Figure B3 the relationship of c, , 0 is 

established. 



cos (0/2) = c/2 B(40) 

R tan 



or 

c = 2 Ro tan 'f cos (0/2) B(41) 

This is the condition £or the camber line to be tangent at 
the exit. 

Using equations B(32) and B(41), 
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2 tan cos < 0/2) 



B<42) 



sin ^ = 

sin < 0/2 ) 

giving, 

cos y- - 2. sin < 0/2) cos < 0/2) 

= 2 C sin < 0/2 ♦ 0/2) ♦ sin <0/2 - 0/2) ] 

= sin 0 B<43) 

Therefore for the camber line tangency condition 

cos ^ = sin 0 B<44) 



Using equation B(39) and B<44) we obtain the following 

2 sin < 0 /2) cos < 0 /2 ) = tan (0/2) 

2 cos 2 ( 0/2) 



and 



2 sin ( 0/2) cos < 0/2) 



2 cos 2 ( 0/2) 



Ro sin 
Ro cos y- 



Ri 



Using the general trignometrical relationship, 

sin 2x = 2 sin x cos x 



we can manipulate the preceeding equation to obtain 

sin 0 [ Rr^ cos y' - R< ] = 2 cos 2 (0/2) 

Ro sin 

Using equation B(43> this becomes 

cos ^ C Ro cos V - Ri ] =2 cos2 < 0/2) 

Using the relationship 

cos2 X = sin2 x cos 2x 

2 cos2 <0/2) = 2 Csin2<0/2) ♦ cos 0] 

= 2 C l* cos 0 ] 

2 



B<45) 



B(46) 
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= 1 ♦ <1- sin2 0) ^ 

Using equation B<43) this becomes^ 

2 cos2 (0/2) = 1 > (1 - co&20'>^ 

= 1 > sin B<47) 

This gives 

Ro cos2 y. - Ri cos y = Ro sin y * Rq ein^y, B<48) 

Writing s^sinO'^ equation <4S) becomes, 

Ro<l - s2) - Rj^ ( 1-32)*^ = Ro 3 ♦ Ro s^ B(49) 

or 

Ro. t <l-s2) - s2 - a] = -d-s2)’^ 

Ri 

(Ro.)2 [ 1 - 2s - 3s2 + 4s 3 ♦ 4s4 ] = l-s2 

(Ri> 

Thus equation B(48) gives the condition 

F<s) = (R< )2 <l-a2) - 1 ♦ a(2 ♦ 3s - 4s2 - 4s3 ) = 0 B<50) 

Ro 

The condition of tangency of the camber line is given 
by the solution of equation B<50) for s between O and 1, 
since s=sin ^ • Using Newton's method of interation the 
root can be found for a given radius ratio. Figure B4 shows 
F(s) in equation B(50) for various radius ratios. Figure B4 
clearly identifies two possible roots. The larger root is 
the correct root^ giving zero blade thickness. The smaller 
root if incorrect, but gives positive values of a and finite 
blade thickness. Emphasis was placed on finding the first 
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root. The value found will give the selected value of sin^. 
Once ^ is founds using equation 6(39) 0 can be calculated 

and c is obtained using equation B<32). 

Using equation B(28) and B<9> 0s and Rs can be 
calculated and the thickness is then given using equation 
B<29) . 

A test is made to assure that the length Ap is 
positive. If Ap is negative this implies a convex rather 
than concave pressure side. This would also imply that the 
inner radius is too small and should be increased to produce 
concave curvature on the pressure side. Ap is found, from 
equations B<19> and B(16). 

The alternate procedure did not require the camber line 
to be tangent at the exit. In this case the value of is 
not solved for but is left as an input parameter. The 
values of 0 and c are then found directly using equation 
B<39) . 

The centers of curvature are given by. 



Xd3 =_Rc. " 


( Rr- 


cos 


( 0/2) 


- Re cos < 0s/ 2) 


) 


cos 


( 0/2) 


Ro 


Ro 






Ro 








yds = Ri 


( Rr^ 


cos 


< 0/2) 


- Re cos < 0S/2) 


) 


sin 


( 0/2) 


Ro 


Ro 






Ro 








XdP = Ec. * 


( Rd 


cos 


< 0p/2) 


- Ec. cos < 0/2) 


> 


cos 


< 0/2) 


Ro 


Ro 




Ro 








ydP = Rl. - 


< Re. 


cos 


( 0p/2) 


- R£i cos < 0/2) 


) 


sin 


< 0/2) 



Rq Ro Ro 

with 
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xd = Rc 

- Rq/Rq 

yd = Ri 

- Ri / R q 

To draw successive blades on a rotor the program must 

rotate the axes about the rotor center after each blade is 

drawn. The following equations were used (see Figure B5> . 

X = X Y sin 0 

cos 0 cos 0 

Y X sin 0 = y 

cos 0 cos 0 

therefore, 

X = X cos 0 + y sin 0 
Y = y cos 0 - X sin 0 

Given the number of blades, the following operations 
(for all three arcs) are performed^ the coordinates are 
calculated as 

^ = xdfds#dP • R cos 0 

y = yd/ds^ds R ain 0 

The arcs are then drawn from 0 = ( (0/2) - (0s, p/2) ) to 

((0/2) *»• (0s, p/2)) 

using 

X = X cos (N) 360 ♦ y sin (N) 360 

2 2 

and 

Y = y cos (N) 360 - x sin (N) 360 

2 2 

where N is the number of blades. Figure B6 and B7 
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illustrates the computer generated shapes and passage for a 



60 blade design. 

B.2 COMPUTER PROGRAM •*ACIRC” 

B.2.1 Description 

A computer program, •“ACIRC** was written to 
construct double circular arc blading using the above 
procedure. The program will calculate and draw blading for 
a initially radial -outflow using Newtons Method to solve for 
a suitable value of 1^, or for an arbitrary input value of y* . 
In either case the method will yield a suction surface which 
is tangent to the outer rim. 

The program is user friendly and is annotated for 
ease of understanding. The program is divided into six 
major parts. 

The first part is the input data section. The 
second part calculates the periodicity, radial solidity, 
spacing, and asks the user if the value of ^ is to be 
calculated or to be specified. 

The third part calculates the first root for 
tangency conditions at the exit. No input is required for 
this section. The fourth part calculates the parameters 
needed for the suction surface to be tangent to the outer 
rim. Usually, if the user specified a value of ^ there is 
the possiblity of the distance from the chord to the 
pressure circle being negative. If this occurs the result 
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is a convex arc. To remedy this, the program asks the user 
to increase the dimension on the inner radius by a 
percentage. The program will use this value to solve for 
the required inner radius to achieve a suction surface that 
is tangent to the outer rim. 

The fifth part prints the output. The sixth part 
draws the blade surfaces. 

B.2.2 Required Inputs 

The following inputs are requested: 

Outer radius 
Inner radius 
Total number of Blades 

Yes or No response for solving for the value of 
Percentage increase in the inner radius - may 
be required, when the value of is specified 
Number of Blades requested to be drawn 

B.2.3 Listed Outputs 

The following outputs are listed: 

Outer radius 
Inner radius 
Number of Blades 
Angle ^ 

Angle 0 

Results of test for tangency 
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Chord length 



Angle a 

Thickness to chord ratio 
Radius of camber line 
Radius of suction line 
Radius of pressure line 
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B.2.4 Program Listing 



Listed hare is the program listing, "ACIRC" for 
the design of double circular arc blades 



10 1 JUNE 19S4*^****^^*^^ 

20 REM******^PROGRhM to ChLCULRTE ?< IiRhU RRC BLRIUNG******* 
30 REM******F0R RN INITIRLLY OUTFLON CRSCRDE****** 

40 REM****^(-*LISING neutons method to solve for the******** 

50 REM****0F RSI********* 

60 D E b 

70 SCRLE -1 j 1 > -1 ’ 1 

80 DISP "SET PLOTTER 7(X) BY 7 < Y ) -<C0NT > " ; 

90 STOP 

95 REM*****IHPIJT BRTR************************ 

100 DISP "ENTER THE OUTER RRDIU3"5 
110 INPUT R0 

120 DISP "ENTER THE INNER RRDIUS"; 

130 INPUT R1 
140 R4=R1/R0 

150 DISP "ENTER THE NUMBER OF BLRDES"; 

160 INPUT 20 
170 L=0 

180 REM********END INPUT DRTR**** 

135 REM*****CRLCULRTIONS********* ************ 

190 S5=360/20 
200 P 1=3. 141592654 
210 S2=2*P1*R0''Z0 
215 R9=C.R0-R1>/S2 
220 S1=<R0-R1)/S2 
230 S0=32/R0 
240 R2=S0*S1 

250 DISP "NRNT PSI CRLCULRTED Y=1>0=N"> 

260 NR IT 1000 

270 INPUT 01 

280 IF 01=1 THEN 340 

290 DISP "INPUT PSI"; 

300 INPUT S3 
310 02=C0SS3 

320 T0=RTN<O2/SQR(1-O2t2) > 

330 GOTO 380 
340 GOSUB 1350 

350 REM*****N0N DETERMINE THE VRLUES FOR PSI>PHIi.RND C***** 
360 S3=RTN<K9/SQR< 1-K9t2>) 

370 T0=2*RTN<R0*SINS3/(R0*COSS3-R1 )) 

380 05=3 1 NS3 
390 O6=CQST0 

400 Cl=R0*SINS3/SIN<T0/2) 

405 REM************************************** 

410 P.EM******TEST FOR TRNGENCY***** 
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420 C 2 = 2 * T a H S 3 * C 0 3 < T 0 ■ 2 > ^ R 0 

430 F.;EM******CCiHCLUriES CRLCULaTION FOR PS I .« PH I - HUH 

440 a 0 = 9 01 + S 3 - T 0 

450 Cti=Cl ''R0 

4S0 R3=C0/ ( 2*S I H < T0/2 ) ) 

470 R4=R1/R0 
4S0 D0=RTH<R3/R4) 

4S5 REri*******'^*'^^'^'^*'^*'********”**-^********* 

490 REM*^*^FOR SUCTION SURFRCE TaHGENT TO OUTER 

500 T2=<2^aO)+T0 

510 R 5 = C 0 ••••■ 2 ^ S I N < T 2 / 2 > > 

520 T 1 = ( ( R 5 * C 1 - C 0 S T 2 / 2 > ) > - '■ R 3 * < 1 - C 0 S 0 T 0 2 ) ) > ) / ( C O • •• ■ 2 ') 

5 3 0 a 1 = R 3 * (1 - C 0 S ( T 0 ■ ■••• 2 > ) + C 0 ^ T 1 .•• ■ 2 
540 a2 = ai-(C0^Tl > 

550 IF a2>0 THEN S40 
560 IF L=1 THEN 620 

570 DISP "NO SOLUTION-RI IS TOO SMaiL"; 

580 nisp " 

590 DISP "H-INCREaSES RI I.E.=25.'; ENTER H"5 
600 INPUT H 
610 L=1 

620 R1=R1+(H/100> 

630 GOTO 200 

640 R6= ( a2/2> + ( C0t2/ < S*a2 ) ) 

650 P4=C0/<2^R6) 

6 6 0 T 3 = 2 * a T N C P 4 / S Q R ( 1 - P 4 1 2 > ) 



fc .70 

680 

690 

695 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

845 

850 

860 

870 

880 

890 

900 



0 3 = R 3 R 0 

O4 = P'5*R0 

0 1 = R 6 ^ Rl 0 

REM*******PRINT STaTEMENTS*^^^^^^^***^^*^ 
PRINT "R0="R0"RI="R1"NO. BLaDES="20 
PRINT "PSI="S3 
PRINT "PHI="T0 

PRINT "CONDITION FOR TRNGENCY IS" 

PRINT " SIN PHI=C08 PS I" 

PRINT " SIN PHI ="05" COS PS I ="06 
PR I NT " CHORD= " C 1 " CHORD TEST= " C2 
PRINT "aLPHH="a0 
PRINT "RRDiaL SOLIDITV=" , R9 
PRINT "THICKNESS TO CHORD RRTIO="iTl 
PRINT "RRDIUS OF CRMBER CIRCLE ="03 
PRINT "RRDIUS OF SUCTION CIRCLE ="04 
PRINT " RRDIUS OF PRESSURE CIRCLE="01 
REM******END PRINT STRTEMENTS***** 



DETERMINED- 



E M * 4 C E N T E R S C R N N 0 W B E 
N 1 =R3- ( < R3*C0S ( T0/2 ) ) - ( R5 
Y1=R4-K CR3^COS(T0^'2) >-(R5 
X2 = R3-K ( R6*C0S ( T3/2 > ) - ( R3 
Y 2 = R 4 - ( ( R 6 * C 0 S ‘s T 3 2 ) ^ ~ ^ F! 3 * C 0 3 k T 0 
N3 = R3 



CO 

CO 

CO 



::''T2/'2) > >*COS(T0 



5-:;T2/2>> 
3 >:T 0 ^' 2 >) 
2 )) 



SINCT0. 

COS(T0. 

SIN<T0 



' 2 ) 
''2 > 
- 2 ) 
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910 Y:5 = R4 

920 REM^^^^DRFiU BLRDE BURFRCES**^* 

930 nrsp "INPUT THE HUMBER OF BLRDE3 TO DRRN"; 

940 NR IT 1000 
950 nisp " 

9fc‘.0 niSP "STRRTING RT H=l"; 

970 INPUT N 

9B0 REM EDIT NEXT LINE FOR PLOT RESOLUTION- 

990 N9 = 25 

1000 FOR N0=1 TO N 
1010 I 1 = < T 0 / 2 > - ( T 2 /■ 2 > 

1020 I 2 = ( T 0 2 ) + 's T 2 ■•••■ 2 > 

1030 I3=a2-Il).--N9 

1040 REM^^^FOR THE SUCTION SIDE^^* 

1050 FOR T4=I1 TO 12 STEP 13 
1 0 6 0 ! 4 = X 1 - ( R 5 ^ C 0 S ( T 4 ) ) 

1070 Y4=Y1+(R5*SIN<T4> ) 

1 0 S 0 ; 9 = X 4 * C 0 S < ( N 0 - 1 > * S 5 > + ( Y 4 * S I N ( < H 0 - 1 > ♦ S 5 ) ) 

1 0 9 0 Y 9 = Y 4 * C 0 S ( ( N O - 1 ) * S 5 > - < X 4 * S I N ( ( N O - 1 > ^ S 5 ) > 

1100 PLOT X9»Y9 
1110 NEXT T4 
1120 PEN 

1130 F! E M * * + * F 0 l: R M B E R L I N E * 

1140 FOR T4=0 TO TO STEP T0/N9 
1150 X 4 = X 3 - ( R 3 * C 0 S ( T 4 > > 

1 1 o 0 Y 4 = 'v' R 3 I N T 4 ) 

1170 X 9 = X 4*C 0 S ( ( N 0 - 1 > * S 5 > + ( Y 4*SI N < < N O - 1 > ^ S 5 > > 

1 ISO Y9=Y4*C0S< (N0-1 > *S5 > - ( X4*S I N < (NO-l )*S5> ) 

1190 PLOT X9»Y9 
1200 NEXT T4 
1210 PEN 

1220 REM^*.**FOR PRESSURE SIDE**** 

1230 l4=T0/2-T3/2 

1240 I5=T0/2+T3/2 

1250 IS=( I5-I4;'/N9 

1260 FOR T4=I4 TO 15 STEP IS 

1 270 X4=X2~ ( R6*C0S ( T4 ) ) 

1 2S0 Y4 = Y2+ R6*S I N ( T4 > > 

1 2 9 0 ! 9 — X 4 * C U 3 ( ' . N y ~ 1 .^ * '3 5 * Y 4 * 3 1 N C ' . N y ~ 1 !^ * 'z- 5 j 

1 3 0 0 Y 9 = Y 4 * C 0 S < ■: N 0 - 1 > * S 5 ) - X 4 * S I N < < N 0 - 1 > * S 5 > 

1305 PLOT X9)Y9 

1310 NEXT T4 
1320 PEN 
1330 NEXT N0 

1340 END 

1350 F:EM*****NENTONS NILL BE USED TO CRLCULRTE^* 
1360 REM*****TRNi:ENCY root BETNEEN 0 £< 90 DEGREE 
1370 REM*****SOLVING FOR THE 1ST ROOT****^^**** 
13S0 T 9=0. 00001 
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1390 K9=0.1 

1400 D E F F N F ( K 9 ) = < 1 - R 4 T 2 * <1 - K 9 f 2 > - K 9 ^ ( 2 + 2 ^ K 9 - 4 ^ 9 1 - 
1410 n E F FHD •: K 9 > = - < E * K 9-2* K 9 ^ R 4 T 2 + 1 2 * K 9 1 2+16* K 9 1 3 > - 
1420 K1=K9 

1 4 3 0 f : 2 = F H F ( K 9 ) / F H D ( K 9 > 

1440 K9=K1-K2 

1450 IF RBS(K2)>fiBSa<9*T9> THEN 1420 
1460 RETURN 



4*K9t3> 



76 



TABLE B1 . List, of Computer and Text Symbols 
for Double Circular Arc Blades 



ACIRC” 


TEXT 


IDENTIFICATION 


AO 


a 


efflux angle 


A1 


As 


suction line to chord length 


A2 


Ap 


pressure line to chord length 


CO 


c** 


c/Rq dimensionless parameter 






for chord length 


Cl 


c 


chord length 


C2 


c 


chord length for tangent test 


11 


- 


Interval ( 0/2 - 0s/2) 


12 


\ 


Interval ( 0/2 * 0s/2) 


13 


- 


Interval (0/2 * 0s/2) ~ (0/2 
- 0S/2) 


K1 


- 


constant value 


K2 


- 


Newtons method, 


K9 


- 


initial guess for Newton^'s 
method 


L 




counter 


N 


N 


Number of blades to draw 


NO 


- 


Number of blades to draw for 
loop 


N9 


~ 


Number of points for interval 
of arc 


01 


- 


Equals 1 for solution of 



for input ^ 
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TABLE Bl. List of Computer and Text Symbols for 
Double Circular Arc Blades (con't) 



;iRC“ 


TEXT 


IDENTIFICATION 


02 


cos ^ 


cos 


03 




output print of Rc 


04 


Ra 


output print of Rs 


05 


sin 


for tangent test 


06 


cos 0 


for tangent test 


PI 


n 


n 


P4 


- 


constant^ cRo/2Rp 


RO 


Ro 


outer radius 


R1 


Ri 


inner radius 


R2 




density^ (Ro " Ri>/Ro 


R3 


Rc /Rq 


camber radius, dimensionless 


R4 


Ri/Ro 


inner radius, dimensionless 


R5 


Rg/Ro 


suction radius, dimensionless 


R6 


Rp/Ro 


pressure radius, dimensionless 


SO 




periodicity, 2 TTRo/Z 


SI 


ffr 


radial solidity 


S2 


s 


spacing 


S3 


f 


blade curvature with respect 
to camber line 


S5 


ys 


periodicity, 360/2 


TO 


0 


camber angle 
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TABLE Bl. List of Computer and Text Symbols for 
Double Circular Arc Blades (con't) 



ACIRC" 


TEXT 


IDENTIFICATION 


T1 


T 


thickness to chord ratio 


T2 


0S 


suction angle 


T3 


0p 


pressure angle 


T4 


d 


loop parameter to draw arcs 


T9 


- 


tolerance on Newtons method 
of integration 


XI 


Xds 


suction center x coordinate 


X2 


XdP 


pressure center x coordinate 


X3 


Xd 


camber center x coordinate 


X4 


- 


variable for x values in 


• 




continuous loop 


X9 


- 


variable for x values in 
rotational loop 


Y1 


Yds 


suction center y coordinate 


Y2 


YdP 


pressure center y coordinate 


Y3 


Yd 


camber center y coordinate 


Y4 


- 


variable for y values in 
continuous loop 


Y9 


- 


variables for y values in 
rotational loop 


ZO 


Z 


number of blades 
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Figur® Bl. Double Circular Arc Blade 
Geowetr ical Description 
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Figure B2. Geometrical Description for Camber Line 

Tangency Condition. 
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Flgura B3. Continued Geoaetrlcal Doocriptlon for Caaber 

Line Tangoncy Condition 
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Figure B4. Function of Resulting fron the Condition 
of Tangency of the Camber Line (equation B(50)). 



S3 




Figur* B5. Scheaiat.ic for Rotation of Axes. 
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Rq s 4 Ri = 3.2 No. Bladaa = 60 

i' - 9.1.72521902 0 » 80.82738092 

Condition of Tangency is 
ain 0 = coa 

ain 0 = 0.159407755 coa » 0.159409430 
chord * 0.983540343 chord teat » 0.983540073 
alpha > 18.34514098 
radial eolldity = 1.909859317 
thickness to chord ratio » 0.194933797 
radius of camber circle ■ 0.758551242 
radius of suction circle » 0.575175292 
radius of pressure circle * 1.462802226 



Figure 86. Computer Generated Shapes 
ioT a 60 Blade Design. 
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Flgur* B7. Co»put«r G«n»rot«d Bxad* 
for a 60 Blade Deelgn. 



APPENDIX C 



WEDGE-ARC BLADING DESIGN 



C.l METHOD 

Because of the limitations seen in the outlet passage 
shapes which could be generated using purely double circular 
arc blades^ symmetrical linear surface segments were added 
in an attempt to achieve a convergent nozzle at a shallow 
angle at the rim exit. 

From the design using only double circular arcs it was 
found that by using the original camber line and the suction 
line to describe the blade^ improved convergency at the exit 
was obtained (see figure Cl). Removing the original 
pressure line would cause the flow to be at a small 
incidence angle at the inlets but this was thought to be 
acceptable . 

By allowing the flow to exit over a straight suction 
side from the exit of the passage^ better control of the 
outlet flow angle was assured. The wedge arc geometry is 
shown in figure C2. 

Referring to figure C2^ the wedge arc blade shape was 
developed to have the following features^ or constraints: 

1. The pressure aide is a circular arc. 

2. The auction side between point H and I is a circular 
arc^ between straight and symmetrical line segments. 
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3. A wedge angle (() ie specified at the leading 
and trailing edges. 

4. Equal leading edge and trailing edge thicknesse ( S) 
are specified. 

5. The radius to the auction surface at the point of 
tangency to the straight segment, (point I) 
intersects the pressure surface of the adjacent blade 
at the trailing edge (point C' > , or downstream of 
point C' p depending on the value of the parameter 

( i'y 

Referring to figure C2 and C3 the following 
relationships are obtained: 



XF = xc - 


S' 


cos 


0 


C(l) 


YF = yc 


S' 


sin 


0 


C(2) 



where S' is the leading edge and trailing edge thickness. 
The value of S' has been made dimensionless by dividing by 
Rq- All other quantities have been made dimensionless in 
the procedures given in appendix B. 





= Xc - 


S/2 cos 0 


C(3) 


yt 


= yc ♦ 


S/2 sin 0 


C<4) 




Xe = 


- 5 


C<5) 




ye = 


Ri 


C(6) 




XL = • 


- s/2 


C<7) 




YL = 


Ri 


C<8) 



Sd 



The wedge angle^ ( is specified. With ( known, the 
length of the straight section i Jl '> from the exit can be 
found. For a symmetric blade, the length (^ ) is the same 
at the inlet and the exit. Referring to figure C2, C3 and 
C4 

^ ' <1 - C(9) 

where Jl - Jl ' when the point I lies on the radius C'D' 

(figure C4) . The slope of the line C'F ( p) is given by, 

p = arc tan ( vr - ) C(IO) 

( xp - Xc^> 

The coordinates of C' are found by solving for the 
previous blade and rotating the coordinates. These values 
are placed in an array until needed. A straight line is 
drawn from the current blade to the previous blade as is 
defined as 

C'F * ( [yF - yc3^ ♦ txF " Xc'3^ ) ^ C(ll) 

Referring to figure C3, the slope of FG is given by, 

A = tot - y’ - ( ] C<12) 

A relationship of X ' f C'F and the line FG is given by, 

X' = C'F cos C'FG 

= C'F cos C oc - - C ♦ p] C(13) 

From equation C(13) the coordinates of I and H are found, 

XI = XF - i cos Ca - ^ - (] C<14) 
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yi = YF - 


sin Ca 


- ^ 


- n 




C<15) 




XH = - 


C 5 + i 


sin 


n 




C(16) 




YH = Ri/Ro * Jl 


cos 


( 




C<17) 


The 


centers of the 


suction 


side 


are given 


by 


either. 




XD' 


= XH ♦ Rz 


cos 


( 




C<18) 




YD' 


= YH ♦ Rz 


sin 


( 




C(19) 


or. 


xd' = XI 


sin [ 


a " 


^ - n 




C<20> 




YD' = YI - 


R^ cos C 


a " 


^ - a 




C(21) 


From 


equations C(17), C(19) 


and 


C(ll) Rz 


can 


be found 


since 


XI - XH 


= R^ (cos 


i - 


sin ^ ) 






yielding. 


Rz 


= XT " 








C<22) 






cos i - 


sin 









To calculate the position of I for a given throat 
dimension <a') at the exit^ further manipulation is 
required. Referring to figure C4, where "a" denotes the 
center of the exit throat, 

Xa = Xc' a' sin a C<23) 

and 

ya = Vo' ~ CO® ^ C(24) 

From equations C(21), C(13) and C(15), 

^2 " < xp - j cos A s 4. jl sin i ) 

(cos i - sin 

= XF ♦ $ ♦ Csin i - cos A 3 C(25) 

Ecos C-sin 
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Referring to figure C4 

[ Rz * a']2 = [ XD'~ xa ] 2 ♦ [ ya - YD'^^ C(26) 

Using equation C(20) and C(21) the above equation becomes, 

[ Rz + a']2 = [xj-Rz sin ” Xa3^ * tya'YI'Rz cos ^3 ^ 
referring to equations C(14), and C<15), 

t Rz ♦ a'] 2 s [xF - J cos A, * A “Xa3 2 

♦ [ ya ~ yr * /^ain ^ - Rz cos ^ J2 

which reduces to 



2Rza' + a' 2 = <xf “ 


Xa> 


2 + 


<Ya - 


yf)2 


* 


C<27) 


•*•2 jl cos 


A 


<Xa 


- Xf) 


* 2 J 


sin ^ (ya 


- yF> 


+2Rz cos 


A 


<Ya 


- yp> 


* 2Rz 


sin A ^^F 


- Xa> . 



Using equation C<25) the following results can be 



obtained 



a'2 s <xF“ Xa>^ ♦^ya" ^ ^ ^ ^ (Xa~XF> 



(cos ( 



♦2 J. sin A (ya 


-yp) 


♦ 2 / xF ♦ S * 


(sin ( - cos A ] 


' (cos 


( - sin ^ ] 


+ 2 / XF ♦ S + 


(sin ( - cos A ] 


^ (cos 


( - sin A. ] 


+ 2 / XF + 8 


•^(sin ( - cos A ] 


1 ( cos 


( - sin A 3 


3a'2 = ( cos 


( - sin “ 



:xF-Xa> 



C<28) 

<XF - Xa>2 + <ya - YF>^ * Jl ^ '\ 

♦ 2/?cosA(xa - xF > ♦ 2 ^sinA <ya-yF)/ 



+ 2 
+ 2 


XF 

xF ♦ 


8 ♦ 
8 + 


^(sin 
J? (sin 


( - COS A 2. 

( - cos A ] 


sin A (xp-x^) 
cos A ( ya-yF) 


+ 2 


xF 


8 ♦ 


(sin 


( - cos A 3 


a' 
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and on combining berms 



2 

Ceos ( - sin ^ la' = C cos ( - sin ^ ] C<29) 

♦ (2 Ceos € - sin ^ 3 cos ^ <Xa - xp) 

♦2 Ceos C - sin A3 sin A (y^ - yp) 

♦ 2 xp sin A(xp - x®> + 2 S sin A(xp -Xq) 

+ 2 Ceos i - sin A 3 sin Ai. (xp - Xa> 

♦ 2 xp cos A<ya ~yF> +25 cosACya ~ YF> 

♦ 2 Ceos C - sinA3 cosA.(ya “ YF^ 

♦ 2 xp a' 2 S a' 2 a' Ceos ( - sin A 3 

Ceos ( - sin 3 a'2 = C cos ( - sinA3 ^2 COO) 

♦ { 2 Ceos € -sin A3 cos A (Xa ” xp) 

♦ 2 Ceos ( - sin A3 sin A (■Ya~ YF^ 

♦ 2 Ceos C - sin A3 sin A<xp- Xa) 

♦ 2 Ceos (. - sin A3 cos A(ya ~YF^ 

+ 2 a' C cos ( - sin A] ) 

2 xp sin A (xp - Xa3 +25 sin A<xp - Xa) 

+ 2 xp cos A<Ya “ YF^ +25 cos A(ya - YF^ 

♦ 2 xp a' *2 S a' 

The above equation is quadratic in the unknown H , and a 
solution can be written explicitly. The value o£ Jl can be 
found using equation CO). The result is as follows. 

If 

B = C 2 Ceos ( - sin A 3 cosACxa ~ xp) COD 
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2 


Ceos 


( - sin A] 


sin A <ya 


- yp> 


♦ 2 


Ceos 


( - sin A2 


sin A < xp 


— Xg ) 


♦ 2 


Ceos 


i - sin A ] 


cos A <Ya 


- yp> 



2 a ' [ cos i - sin A] J / Jf '^cos i 



Let C = 



then our 



2 xp sin A <xp -Xq) + 2 S sin ^<xp - Xg) 

2 xp cos tK^Ya ~YF^ ♦25 cos A <ya “ yp) 

2 xp a' ♦ 2 5 a'}/ ^ [cos ( - sin 3 

quadratic equation is. 



O 



J^'^Ccos ( - sin 

<1 - r ) 



_] 2B J?'^[cos 

♦ C "^(cos { -sin 3 



-sin ] 



and is of the form, 

O = 



Z 

( 1 - (' ) 



2 Bd-(') ^ C 



so that 

<!-[') = -2B/2 ± V(2B)^ - 4C 

2 

= - B + V< B C) 

The following possibilities occur: 

1. If the term under the square root is negative 
solution does not exist. 

2. If the term under the square root is positive 

the value of (1 - (') is between 0 and 1 then 

solution existsi the solution is 



( 1 - (1 



-B ♦ KB^ - C) 



sinA 3 
C<32) 



» C(33) 
C<34) 

C(35) 



C<36) 
then a 
and 
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3. If the value of <1 - Cl) obtained in step <2) is 
outside the range of 0 to 1 then the solution is, 

< 1 - (1) = -B - C) 

4. If the value of (1 -Cl) obtained in step (3) is 
outside the range of 0 to 1 then no solution exists. 

Once C' has been determined and all other values have 
been found, the design is complete and the blades can be 
drawn . 

Referring to figure C2, to draw the radius on the 
leading edge the points 



XL ♦ 


S/2 cos 


6 


C<37> 


YL - 


S/2 sin 


e 


C<38> 



are plotted over the interval of 9 between O and 180 
degrees. The radius on the trailing edge is given by. 



XT - 


S/2 cos 


d 


C<39> 


YT 


S/2 sin 


e 


C(40> 



plotted over the interval of © between 0 and 0-^180 . 

To draw successive blades on the rotor the procedure is 
to rotate the axes after each blade is drawn. The procedure 
is described by equation B<55) and B(56) . 
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C.2 COMPUTER PROGRAM ••AWEDG'* 



C.2.1 Description 

The wedge arc blade design program, "AWEDG** 
calculates and draws wedge arc blading for an initially 
radial -outf low rotor. The difference between '*AWEDG** and 
“ACIRC** is not only the specified input parameters; ''AWEDG** 
calculates the straight section needed on the auction side 
to constrict the flow to a know dimension at the trailing 
edge tip of the adjacent blade. The program, for simplifi- 
cation, is divided into seven parts and each part will be 
described . 

The first part (lines 1 to 460) contains the 
values of all specified parameters. On running, it also 
allows the user to specify the angles ^ and 0 to determine 
early in the program what the value of a will be. If the 
user wishes to change the value of a he can do so. Chord 
length, camber radius and a are calculated in this section. 

The second part (lines 470 to 640) prints some 
of the output parameters and calculates the center 
coordinates of the camber line. 

In order to calculate the outlet passage size 
the preceding passage must be analyzed. The program does 
this in part three (line 650 to 800) by calculating the 
position of the trailing edge tip of the preceding pressure 
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line. Point C' on figure C3 is calculated by rotation of 
axes . 

In part four (lines 810 to 1020) the position of 
points C, F, T, E, and L, shown in figure C3, are 
calculated . 

Part five (beginning on line 1030) starts the 
plotting routine- The number of the blades requested to be 
drawn are specified- The pressure side radius is calculated 
and is drawn. Calculations for the straight line segments 
are performed (lines 1340 to 1570). Based on the determined 
value of the straight segment the coordinates of point I are 
dervied (lines 1630 to 2090). 

In part six the radius of the suction side is 
drawn. Plotting for the leading and trailing edge radius 
are also accomplished in part six. 

Part seven (2100 to 2140) concludes the 
procedure by printing the value of the throat diameter, 
length of the straight line segmnent and radii of pressure 
and suction circle. 

C.2.2 Required Inputs 

The following inputs are required: 

Outer radius 
Inner radius 
Total number of blades 
angle ^ 
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angle 0 (must be consistent with -see 
discussion) 

Diameter of the circle at the passage exit 
Leading edge and Trailing Edge thickness 
Leading edge and Trailing Edge wedge angle 
Number of blades requested to be drawn 

C.2.3 Listed Outputs 

The following quantities are output: 

Outer radius 
Inner radius 
Number of blades 
Angle ^ 

Angle 0 
Angle a 

Radial solidity 
Chord length 

Thickness to chord ratio 

Leading edge and trailing edge thickness 

Leading edge and trailing edge wedge angle 

Throat width at exit 

Radius of the pressure surface 

Radius of the suction surface 

Length of the straight section on the 

suction side. 
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C.2.4 Proqraa Listing 



Liated here is the program "AWEDG" for the 
design of wedge shaped are blades. Table C2 identifies the 
symbols used in the program with the defined parameters. 
Figure C6 and C7 illustrates computer generated shapes for 
an outer radius of 4 inches, an inner radius of 3.2 inches 
and 60 blades. The value of ^ and 0 are specified 
parameters to limit the efflux angle a to between 8 and 9 
degrees . 



10 R E M ^ I’i E D G 1 fi U G U S T 1 9 y 4 ^ 

20 REN**^^PROGRflM TO CflLCULflTE S. DRRH ARC BIRDIHG^** 

30 REM^^^RRCS RHD STRRIGHT BROKING FORI RN*** 

35 REM***IHITIRLLY OUTFLOW CRSCRDE^^-^^^*^ 

40 DEG 

50 DIM GC 100 3>HC 10U]>UC 100]j VC 100] 

G 0 S C R L E - 1 j 1 j - 1 j 1 

70 DISP "SET PLOTTER 7(X> BY 7 < Y> -< CONT > " j 
80 STOP 

100 DISP "ENTER THE OUTER RRDIUS"; 

110 INPUT RO 

120 DISP "ENTER THE INNER RRDIUS"; 

130 INPUT R1 
140 R4=R1/R0 

150 DISP "ENTER THE NUMBER OF BLRDES"; 

160 INPUT 20 
170 L=0 
130 S5=360/20 
190 PI =3. 141592654 
200 S2=2*Pl+R0/20 
210 R9 = (R0-Rn''S2 
220 S1 = '::R0-R1>/S2 
230 S0=S2/R0 
240 R2=S0*S1 

250 DISP "ENTER PSI-< BTWN LE & TE 0-360/'2"; 

260 INPUT S3 

270 DISP "ENTER PHI -PRESSURE RNGLE/CRMBER< " ; 

230 INPUT T0 
290 R0=9O+S3-T0 

300 PRINT "RLHPR IS DESIRED BWTN 8-9 DEGREES. THE VRLUE 
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7tj REM: 



310 
320 
330 
340 
350 
360 
370 
330 
3 9 0 
400 
410 
420 
430 
440 
450 
460 
470 
430 
490 
500 
510 
520 
530 
540 
550 
560 
570 
530 
590 
600 
610 
620 
630 
640 
650 
660 
670 
630 
690 
700 
710 
720 
730 
740 
750 
760 
770 
730 
790 
3001 



PRINT "OF PSI!1:PHI CHOSEN WILL YIELD RN" 
PRINT "RLPHR="R0" DEGREES" 

PRINT "PSI="S3"RND PHI="T0 

DISP "NISH TO CHRNGE P3IS:PHI Y=1jN = 0"; 

INPUT 01 

IF 01=1 THEN 250 

DISP "ENTER THE LE S. TE THICKNESS"; 
INPUT Q1 
M0 = Q1 •••■RO 

DISP "ENTER LE S-; TE HEDGE RNGLE"j 
INPUT Ml 

C 1 = R 0 ^ S I N S 3 / S I N ( T 0 / 2 > 

C0=C1--R0 

R 3 = C 0 •••■■ ( 2 S I N ( T 0 ■ 2 > > 

R4=R1/R0 
U 3 ~ Fl 3 ^ F! y 



REM*******F'RINT STRTENEflTS*^^***^^**- 



" R0= " RO " R I = " R 1 " NO . BLRDES= " 20 
"PSI="S3 
"PHI="TO 
"CH0RD="C1 
"RLPHfi="RO 

"RRDIRL SuLIDITY="R9 
"LE S: TE THICKNESS = "Q1 
"LE & TE HEDGE RNGLE="M1 
ENTER THE THRORT LENGTH R*"; 
INPUT R5 

REM^^**^^END PRINT STRTEMENTS#*^^* 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

DISP 



REM*:^**CENTERS CRN NON BE DETERMI NED^^** 

J--!3 = R3 
Y3 = R4 

REM**^^DRRH BLRDE SURFRCES^*** 

REM***FIRST CRLCULRTE COORDINRTE OF C’ ON BLRDE NO, 
N9=16 
NS = 7 
Nl=20 
1 = 1 

FOR T4=0 TO TO STEP T0.-N9 
K 4 = :■< 3 - ( R 3 ^ C 0 S ( T 4 ') 

Y4=Y3+(R3*SIN(T4>> 

9 = Y 4 * C 0 S < ( t U - 1 > * S 5 ■;• + ( Y 4 ^ S; 1 1 1 ( ( N 1 - 1 ) * S 5 > 

Y 9 = Y 4 * C 0 S ( ( til - 1 > * S 5 > - < X 4 * S I N < ( N 1 - 1 ) * S 5 ) > 

PLOT X9>Y9 
UC I ]=X9 
V[ I ]=Y9 
1 = 1 + 1 
NEXT T4 
PEN 
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S05 REM***F ?< C ON PRE'dbURE LIHE*'^* 



810 


REM** 


*SEC 


OND C 


ALCULATED 


THE COORDINATES 


S£0 


N0=1 










830 


1 = 1 










840 


FOR T4=0 


TO TO 


STEP T0/ 


N9 


8-jO 


X4 = X3 


-(R3 


*i.:08's 


T4) ) 




880 


Y4 = Y3 


+ K F-! o 


*SIN( 


T4:) ) 




870 


X9 = X4 


^COS 


< (N0- 


1 >*S5)+( Y 


4*SIN< (N0-1 >*S5) ) 


880 


' i' 9 = Y 4 


^cos 


(NO- 


1 >*S5>-(X 


,4*SIN( (N0-1 )*S5) ) 



OF* 



39y 
900 
910 
920 
930 
940 
950 
960 
970 
900 
990 
1 000 
1010 
1 020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 



GC I ] = X9 
HE I ] = V9 
1 = 1 + 1 
NEXT T4 
M2-GC N9+1 ] 

M3 = r‘l2-M0*COST0 
h 4 = N 2 - ( M 0 / 2 ) * C 0 3 T 0 
L2 = HCN9 + 1 ] 

L3=L2+M0*SINT0 

L4=L2+<M0.'-2-'*SINT0 

M5=-M0 

M6=-(M0/2) 

L5=R1/R0 
L 6 = R 1 R 0 

REr'1*****PL0TTING SECTION******** 
DISP "INPUT THE NUMBER OF BLADES 
INPUT N 

REM EDIT NEXT LINE FOR PLOT 

N9 = 16 

FOR N0=1 TO N 
REM****FOR PRESSURE 
1 = 1 

FOR T4=0 TO TO STEP 
X 4 = 3 - ( R 3 * C 0 S ( T 4 ) ) 

V4=Y3+<R 
X9=X4*C0 
Y 9 = Y 4 * C U 
GC I ]=X9 



TO DRAM"; 
RESOLUTION- 



LINE* 



T0/N9 



!■ 



SIN':;T4) > 

C N0- 1 > *S5 > + ( Y4*S I N ( ( N0- 1 > *S5 > > 
( N 0 - 1 > * S 5 ) - ( X 4 * S I N ( < N 0 - 1 > * S 5 > ) 



1170 HE I ]=Y9 
1180 1=1+1 
1190 PLOT X9>Y9 
1200 NEXT T4 
1210 PEN 

1220 REM***SUCTION SIDE LINE***** 

1230 REM***VALUE OF L SOLVED FOR MILL BE GENERATED*** 
1240 REM***TO GIVE THE DESIRED THROAT AT EXIT******* 
1250 E1=L3-VE N9+1 ] 

1260 E2=M3-UEN9+1 ] 

1270 L0 = ATN«:-E1/E2> 

1280 C3=SQR(Elt2+E2t2) 
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1290 D1=RQ-S3-M1 
1 3 0 0 B 9 = C 3 * C 0 S < R 0 - S 3 + L 0 - M 1 > 

1310 REM^^^CRLCULRTIUN OF BLRDE SUCTION SUF'.FRCE:^^^ 
1315 REt'1* + *FOR R GIVEN THRORT SIZE RT EXIT^:^* 

1320 REM^+*R4 IS THE RRHIUS RT THE THRORT R5/2^^^ 

1 3 3 0 R 4 = R 5 / i 2 * R 0 ) 

1 3 4 0 X t. =U [ N 9 + 1 ] + R 4 ^ S I H R 0 
1 3 5 0 V 5 = V C N 9 + 1 ] - R 4 * C 0 S R 0 

1360 2 5 = < S I H D 1 * < V 6 - L 3 > + C 0 SHU ( K 6 - 1'1 3 > ) * < C 0 S f1 1 - S I N B 1 > 

1 3 7 0 Z 6 = < S I H D 1 * ( t'1 3 - X 6 > + C- 0 S D 1 * ( V 6 - L 3 > - R 4 ) * ( S I H h 1 - C 0 S D 1 
1 380 0 1 = ( Z5 + Z6 ') / ( C3* < COSM 1 -S I ND 1 ) > 

1390 Z 7 = ( 2 * ( M 3 + f'1 0 ) * ( S I N D 1 ^ ( t'1 3 - X 6 + C 0 S D 1 * ( V 6 - L 3 ) - R 4 ) ) 
1400 Z8= < ( t'l3-X6 > t2+ < Y6-L3 ) t2-R4 12 ) * COSM 1 -S I ND 1 > 

1410 Q2= 'I Z7 + Z3 > / •: C5t2* ( COSM 1 -S I ND 1 ) > 

1420 Q3=Qlt2-Q'2 

1430 "if Q3<0 them 1560 

1440 Q 4 = - Q 1 + S I j R ( 0 3 > 

1450 IF 0<Q4<1 THEN 1470 

1460 IF Q4<0 OR Q4>1 THEN 1500 

1470 R E M * F I R S T S U L U T I U N F 0 R E = '•. 1 ~ L ! 4 

1480 T9=l-Q4 

1490 GOTO 1550 

1500 Q5=-Q1-SQR<Q3) 

1510 IF 0<Q5<1 THEN 1530 

1520 IF Q5<0 OR 05 >1 THEM 156u 

1530 REM^^^SECOMD SOLUTIOM FOR E=a-Q5>*** 

1540 T9=l-Q5 
1550 GOTO 1570 

1560 PRINT "MO SOLUTION EXISTS" 

1570 B8=B9^a-T9) 

1580 M7=M3-B8^C0SD1 
1590 L7=L3-B8*SIMD1 
1 6 0 0 M 8 = - < M 0 + B 8 * S I M M 1 > 

1610 L 8 = R 1 R 0 + B 8 * C 0 S M 1 

1620 REM**^RRDIUS OF SUCTION SIDE*** 

1 6 3 0 R 7 = < M 7 - M 8 > / ( C 0 S M 1 - S I N D 1 > 

1640 REM***CENTERS OF SUCTION RRDIUS*** 

1650 X5=M8+R7*C0SM1 
1660 V5=LS+R7*SINM1 

1670 REM***PLOT ROUTINE TO DRRN SUCTION SIDE*** 

1680 REM***RLOT RRDIUS ON LE TIP*** 

1690 FOR T4=0 TO ISO STEP ISO.- NS 
1700 X 4 = M 6 + ( M 0 ••••■ 2 ') * C 0 S T 4 
1710 Y4 = L6- < M0.'-2 ) *S I HT4 

1720 X 9 = X 4 * C 0 S ( ( N O - 1 > * S 5 > + ( Y 4 * S I N C ( N O - 1 ) * S 5 > > 

1730 Y9=Y4*C0S < ( NO- 1 ) *S5 > - < X4*S I N ( ■: NO- 1> *S5 ) > 

1740 PLOT X9>Y9 
1750 NEXT T4 
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1 7b0 
1770 
1730 
1790 
1 300 
1 3 1 0 
1320 
1330 
1340 
1350 
1330 
1370 
1330 
1390 



K4=t'18 
V4 = L3 

!• < 9 = X 4 ^ C 0 3 ( ( H 0 - 1 ) * 3 5 ) + < Y 4^31 H ( ( H 0 - 
V 9 = V 4 * C 0 3 ( ( H 0 - 1 ) ^ 3 5 ) - < X 4*31 H ( ( H 0 - 
PLOT X9jY9 

REM***PLOT 3UCTI0H RflEiIUS*** 

I l=-r'll 

I2=2*M1+T0 

I 3= ( 12-11 ) /H9 

FOR T4=I1 TO 12 3TEP 13 

X4 = X5-R7*i:03T4 

Y4=Y5+R7*SINT4 



1) 
1 .j 



• 35 ) > 



X9=X4*C03( (NO-1 )* 
Y9=Y4*C0S( (M0-1 



)+<Y4^SIH< (H0- 
H0- 



35)-(X4*SIN-; 



1 ) 
1 ) 



) ) 
^ 'i 



1900 


PLOT X9jY9 


1910 


NEXT T4 


1920 


X4=M3 


1930 


Y4=L3 


1940 


X 9 = X 4 * C 0 3 < ( N 0 - 1 ) ^ 


1950 


Y9=Y4*C03( (NO-l )■! 


I960 


PLOT X9.-Y9 


1970 


REM***PLOT RfiDIIJ;; 


1930 


1 1=T0 


1990 


I2=T0+130 


2000 


I3=< I2-I 1 )/H3 


2010 


FOR T4=I1 TO 12 : 


2020 


X4=M4-(M0/2)*CO3' 


2030 


Y4=L4+<M0/2)#SIH' 


2040 


X 9 = y 4 * c Q 9 ( ( H y - 1 > 


2050 


Y9=Y4*C03( (M0-1 ) 


2030 


PLOT X9»Y9 


2070 


NEXT T4 


2030 


PEN 


2090 


NEXT NO 


2 1 0 0 


PRINT "H*="fi5 


2110 


PRINT "L="BS 


2120 


PRINT "RRDILI3 OF 


2130 


PRINT "RADIUS OF 


2140 


END 



.5> + 'Yt’4*3IH( 
;5>-(X4*3IH( 



Hy- 

:H0- 



1 ) 

■ 1 ) 



*35)) 



OH THE TE TIP- 



3TEP 13 



1 ) 
1 ) 



35) ) 



PRE33URE CIRCLE 
3UCTI0H CIRCLE=' 



= 03 
'•;r7*R0) 
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TABLE Cl 



List of Computer and Text Symbols 
for Wedge-Arc Blade. 



"AWEDG” 

AO 

A4 

A5 

BB 

B9 

CO 

Cl 
D1 
El 
E2 
G CI3 

Htl] 

I 

II 
12 
13 



TEXT IDENTIFICATION 



(X 



2*a' 



r 



c" 

c 



efflux angle 

throat dimension (radius of 
outlet circle) 

throat diameter (diameter of 
exit circle) 

length of straight section 

length of straight section to 
be tangent to outlet circle 

c/Rq dimensionless, 
chord length 

chord length 

constant angle 

constant value 

constant value 

array for current pressure 
side X coordinates 

array for current pressure 
side y coordinates 

counter 

Interval variable 
Interval variable 
Interval variable 
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TABLE Cl 



List of Computer and Text Symbols for 
Wedge-Arc Blade (con't). 



AWEDGE" 


TEXT 


IDENTIFICATION 


L2 


yc 


final value on present 
pressure line^ y coordinate 


L3 


yF 


final value on present 
pressure line plus blade 
thickness^ y coordinate 


L4 


yt 


center point of trailing 
edge, y coordinate 


L5 


yE 


Leading edge initial suction 
side, y coordinate 


L6 


yi 


point of tangency for 
suction side, y coordinate 


L7 


yi 


point of tangency for 
suction side, y coordinate 


LS 


yH 


initial point for suction 
side curve, y coordinate 


MO 


s 


Leading and trailing edge 
thickness 


Ml 


( 


Leading and trailing edge 
wedge angle 


M2 


Xc 


final value on present 
pressure line, x coordinate 


M3 


XF 


final value on present 
pressure line plus blade 
thickness, x coordinate 


M4 


xt 


center point of trailing 
edge, x coordinate 


M5 


X£ 


Leading edge initial suction 
side, X coordinate 
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TABLE Cl 



List of Computer and Text Symbols for 
Wedge-Arc Blades (con't). 



AWEDGE“ 


TEXT 


IDENTIFICATION 


M 6 


XL 


center point of leading edge 


M7 


XI 


point of tangency for suction 
side, X coordinate 


MS 


XH 


initial point for suction 
side curve, x coordinate 


N 


N 


Number of blades to draw 


NO 


- 


Number of blades to draw for 
loop 


N9 


- 


Number of points for interval 
of arc 


01 


- 


Equals 1 to change ^ and 0 , 
0 for no 


03 


Rc 


output print of R^ 


04 


Rs 


output print of R 3 


PI 


n 


n 


Q 1 


- 


constant 


02 


- 


constant 


03 




constant, value under square 
root term for solving for 
(1 


04 


- 


same as Q3 


05 


- 


same as Q3 


RO 


Rq 


outer radius 
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TABLE Cl 



List of Computer and Text Symbols for 
Wedge-Arc Blades (con't). 



AWEDGE" 


TEXT 


R1 


Ri 


R2 


e 


R3 


Rc/Ro 


R4 


Ri / Rq 


R6 


Rp / Rq 


SO 


s 


SI 


O'r 


S2 


s/Ro 


S3 


f 


S5 




TO 


0 


T4 


- 


T9 


- 


X2 


XdP 


X3 


Xd 



X4 

X9 



IDENTIFICATION 
inner radius 

radial chord, (Rq - Ri>/Ro 
camber line radius, 
dimensionless 

inner radius, dimensionless 
pressure radius, dimensionless 
peripheral spacing, 2 ttRo/Z 
radial solidity 
dimensionless peripheral 
spacing 

blade curvature with respect 
to center of rotor 
spacing angle, 360/Z 

camber angle 

loop parameter to draw arcs 
value of 1 - ( ' 

pressure surface center, x 
coordinate 

camber line center, x 
coordinate 

variable for x values in 
continuous loop 

variable for x values in 
rotational loop 
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TABLE Cl 



List of Computer and Text Symbols for 
Wedge-Arc Blades (con^t). 



AWEDGE" 


TEXT 


IDENTIFICATION 




Y2 


YdP 


pressure center 


y coordinate 


Y3 


Yd 


camber center y 


coordinate 


Y4 


- 


variable for y values in 
continuous loop 


Y9 


- 


variables for y 
rotational loop 


values in 


20 


Z 


number of blades 
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Figure Cl. Relaxed Double Circular Arc Passage 
with Pressure Line Removed. 
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Figure C2. Wedge-Arc Shaped Blading Geoaetry. 
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\ 




Figure C3. Enlarged View of Trailing Edge. 
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Not«: C'D'la not noraol to tha proaaure aurfaca. 



Figure C4 . Geometrical Relationahip Between Throat 
Dimension and Suction Surface Radius 




Rq =4 Ri = 3.2 No. of Blades = 60 
= 9.17 0 = 90.68 

chord = 0.896198294 
a 1 pha = 8 . 49 

radial solidity = 1.909859317 
L.E. St T.E. thickness = 0.02 
L.E. St T.E. wedge angle = 5 
A* = 0.05225 
L = 0.037306915 

radius of pressure circle = 0.629980625 
radius of suction circle = 0.476827358 





Figure C5 . 



Computer Generated Shapes for a Cascade of Wedge 
Arc Shaped Baldes with 60 Blades. 
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Flguro C6. Conputer Generated Blade 
for Wedge Arc Shaped Blade. 
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